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Abstract 

Quantum mechanical real-time tunneling through general scattering potentials is studied in the 
semiclassical limit. It is shown that the exact path integral of the real-time propagator is dominated 
in the long time sector by quasi-stationary fluctuations associated with caustics. This leads to an 
extended semiclassical propagation scheme for wave packet dynamics which accurately describes 
deep tunneling through static and, for the first time, driven barrier potentials. 
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I. INTRODUCTION 



Tunneling through a potential barrier is one of the most fascinating aspects of quantum 
mechanics. In recent years a particular challenge has been to understand tunneling in 
complex systems using semiclassical methods. However, any simple description is seemingly 
hampered by the fact that a quantum mechanical object running towards a barrier with 
a typical energy smaller than the barrier height may penetrate it even though all classical 
trajectories with such energies are reflected completely. Thus, for static observables like 
e.g. tunnel splittings one works in the energy domain and calculates the energy dependent 
Green's function semiclassically by switching from a real-time orbit outside the barrier to an 
orbit in imaginary time, i.e. with imaginary momentum, under the barrier. This technique 
can be traced back to the "old" WKB approximation HJ and meanwhile has been successfully 
extended to extract e.g. tunnel splittings also for systems with classically chaotic dynamics 
0. The crucial question is then: Can semiclassical tunneling also be described in the real- 
time domain? This issue has turned into a fundamental challenge for our understanding 
of semiclassics in general and systems with explicit time dependence as e.g. in the context 
of driven tunneling and control of tunneling || or tunneling in the presence of chaos |3J] in 
particular. 

The probability amplitude for a particle initially at qi to be at qf after time t is given by 
Feynman's path integral representation of the propagator as [[5] 

G{qf, qi ,t) = (q f \e- lHt/h \ qi ) = J V[q)e iS ^ h . (1) 

The integral sums over all paths running from q^ to qf in time t where each contribution is 
weighted according to its action S[q] = jjjf dq\jlm\E — V(q)] — Et with potential V(q) and 
energy E = E(qf, qi, t). Due to the oscillating integrand in ([I]) tunneling appears to be the 
result of a complex interference pattern. In the semiclassical limit the path sum is dominated 
by the contributions of the stationary paths 55^] = obeying Newton's equation of motion 
and small fluctuations around them. In the last decade efficient semiclassical propagation 
schemes based on Gaussian wave packets have been developed, certainly the most powerful 
known as the Hermann-Kluk propagator (HK) ||. However, the inclusion of deep tunneling 
(may be even in presence of external driving) has not been satisfactory yet. In fact, it 
was found that classically allowed real-time trajectories running over the barrier are not 
sufficient to capture strong tunneling |7|, §]. Various extensions have thus been attempted. 
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By propagating a large number of initial wave packets tunnel splittings in a double well 
potential could be extracted, however, in single wave packet motion tunneling effects were 
absent pi. As for barrier penetration no low energy stationary orbit exists along the real 
axis, it is tempting to think that one may find one in the complex plane ||. An individual 
path circumventing-in real-time-the barrier region in a complex coordinate plane is assumed, 
but its existence is still elusive. In [ID| orbits run along the real axis, but the semiclassical 
propagation over the time interval t is time sliced into steps over intermediate intervals. This 
spawning of orbits turns out to be numerically extremely expensive already for two slices 
and improves tunneling amplitudes only for energies not too far below the barrier top. 

In the sequel we re-examine the semiclassical barrier penetration through scattering po- 
tentials starting from the exact expression ([!]). Our idea is this: While in the energy domain 
tunneling is described within a complex time plane, here, we study classical mechanics for 
complex energies. Our analysis reveals how tunneling is encoded in the quantum propaga- 
tor ([!]) in terms of real-time orbits. In particular, it turns out that an individual complex 
"tunneling path" does not exist. Based on these results we extend the conventional HK to 
semiclassical wave packet dynamics in the deep tunneling regime not only through static 
potentials, but, for the first time, report also on accurate results for driven tunneling. The 
approach is shown to be very efficient for one-dimensional systems and may thus serve as a 
promising starting point for higher dimensional studies. 



II. COMPLEX MECHANICS AND SEMICLASSICAL APPROXIMATION 

We consider the motion of a particle of mass m = 1 in a general one-dimensional, sym- 
metric barrier potential V(q) where the barrier top is located at q = with V(0) = Vq. V(q) 
is assumed to be a smooth and analytic function of q that can be approximated around q = 
by an inverted harmonic oscillator and for large \q\ falls off as V(q) — > Vo/[q/l} 2k ,k > 2, 
integer, with a typical barrier length scale I. A sufficiently high barrier is taken for granted. 
Typical examples include V^{q) = V /[l + (q/l) 2 ] k , but our results also apply to the Eckart 
barrier Vq/ cosh(g/Z) 2 and the Gaussian barrier Vq exp(— q 2 /l 2 ). 

Now, think of a wave packet ip(qi, 0) localized to the far right (g, > 0) which is propagated 
towards the barrier according to ip(qf, t) — J dqiG(qf, q^ t)il){q^ 0). We are interested in that 
portion of the packet that after time t arrives on the far left (g/ = — g, < 0). All real orbits 
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connecting the two asymptotic regions run over the barrier (E > Vq) and as, for fixed end- 
points, t becomes large they spend most of their time in the parabolic range around q = 0. 
There, the marginal stability of trajectories causes the semiclassical G(—qi,qi,t) to die out 
exponentially in contrast to exact results |7| . Classical orbits with E < Vq coming from the 
far right or left reach the right or left flank of the barrier at turning points (TPs) go an d 
—go, respectively; the long time properties of the path integral ([I]) are therefore governed 
by the dynamics in the "forbidden" range between the TPs. Mathematically, for E < Vq no 
real stationary phase point to (H) obeying the proper boundary conditions exists in function 
space. The usual procedure is then an analytic continuation meaning here to extend classical 
mechanics to the complex coordinate plane. 

Newton's equation of motion, q + V'(q) = 0, where q = dq/dt and V = dV/dq, translates 
for complex q = x + iy into 

x + r x = 0, y + j x = 0. (2) 

Here, V(q) = r(x,y) +ij(x,y) and the subscript x [y] denotes the partial derivative with 
respect to x [y]. We further exploited that for analytic functions V(q) Cauchy's relations 
r x = i y and r y = —j x apply. From (§) one simply finds that the total energy E = e rc + ie im 
and its real and imaginary parts 

e re = {x 2 -y 2 )/2 + r(x,y) and e im = xy + j(x, y) , (3) 

respectively, are constants of motion. How does the corresponding classical mechanics look 
like? For low energies the TPs go, —go lie in the range where V(q) can be approximated by 
its asymptotic behavior. Hence, we consider paths starting from large g^ = x-i > along 
the real axis with complex momentum p$ = q^ = Xi + iiji. Typical trajectories are depicted 
in fig. [I]. While basically three kinds of orbits can be distinguished, the common behavior 
is that as the barrier vanishes asymptotically, for large distances from the top the classical 
motion tends to be a free motion. If we represent trajectories in the form q{t) = i?(t)e 1 ^, 
for very large R they run close to straight lines with constant 4>{qi,Pi) depending merely on 
the initial phase space variables. Let us now discuss the types of orbits in detail. 

First, we look at class (a) as it is the only one where orbits cross the line x = (<ft = 7r/2) 
to reach the other side of the barrier. For such orbits we need initially Xi < (otherwise they 
would run away from the barrier anyway) and for the discussion assume iji > implying 
e im = i-iVi < 0, see Eq. (|). We further note the asymptotic form of V(q) = r(g) + ij(q) 
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using polar coordinates: 

(■d a\ cos(2/c0) sin(2A;0) 

One sees immediately that a successful crossing of the dividing surface must happen with 
y > (and i; < of course). Namely, at = tt/2 the imaginary part j vanishes so that e; m = 
xy = iif/i < meaning y > 0. For the required energy one derives |e im | > r(xi,0) = V(xi). 
Starting, however, from the imaginary axis with x < and y > 0, i.e. in the direction of 
decreasing |V(g)|, always generates an orbit reaching the asymptotic left side of the barrier 
far from the real axis (<f) < tt). We conclude that a simple "tunneling path" connecting 
the asymptotic segments of the real axis on either side of the barrier via a tour through 
the complex plane does not exist. This is in sharp contrast to tunneling for fixed energy. 
There, the energy dependent Greens function K(qj, qi, E) exhibits stationary phase points 
in imaginary time corresponding to classical paths running with energy E in the inverted 
potential through the barrier range from qi to qj. Here, for real-time tunneling a stationary 
phase path to the quantum propagator G(qf,qi,t) even with complex energy cannot be 
found. This important result may also reflect the quite different roles "energy" and "time" 
play in quantum mechanics. 

For our analysis the consequences are two-fold: on the one hand complex trajectories 
in class (a) do not play any role for a semiclassical approximation to G(qj,qi,t), and on 
the other hand the path integral in ([!]) is in the low energy sector completely determined 
by fluctuations. To find its dominant contributions thus means to detect the dominant 
fluctuations; these are points in function space which lie close to orbits with 5S[q] = 
and also obey the proper boundary conditions. Accordingly, we consider the remaining two 
classes of paths. 

The second class (b) contains orbits with small but non- vanishing energies < |ei m | < 
r(xi, 0) which may exhibit TPs in the complex plane and always live on the same side of the 
barrier. Hence they are not relevant either. In the third class (c) trajectories have real total 
energy E, i.e. e im = 0, but start with purely imaginary momenta x« = and, as assumed, 
small e re = —yf/2 +r(xi,0). These orbits display crucial features as we will explain in the 
following. For that purpose we focus on the limit e re = and follow paths with X{ > 0, yi > 0. 
Writing asymptotically q(t) = R(t)e lr ^ one obtains 

<P(x t ,y l )=<f ) t = 7i/{2(k + l)]. (5) 



Hence, after a transient period of time all those orbits run along the same line in the complex 
plane independent of their starting points x,. And since they carry the same energy, they 
are also focused in phase space so that the line C defines a caustic. Were the trajectories 
optical rays, = 0+ would be a burning line. Due to symmetry the same holds true for the 
complementary line — 0+ and the lines 7r ± 0+ on the other side of the barrier. Typically, a 
caustic is associated with unstable orbits and fluctuations connecting them which renders a 
simple Gaussian semiclassics insufficient [[|. To verify this scenario here, we consider small 
deviations Sq = Sx + iSy around a certain orbit q(xf,t). By linearizing the equations of 
motion (||) one gains (Sx, 5y) T =M (5x, 5y) T where M is the stability matrix evaluated 
along q(t). Along = 0+ its diagonal elements are — f xx > and the off-diagonal elements 
vanish j xx = 0. Accordingly, all trajectories merging along the burning lines are unstable. 
Small deviations in phase space can lead from an orbit q(xf,t) to another one q(x'f,t) and 
even allow for a turn from positive to negative momentum to run along the q{x\] i)-orbit 
back towards the real axis. As asymptotically paths creep along + , jumps from very small 
positive to negative momenta require only tiny fluctuations. The reversed orbit crosses the 
real axis at x\ and approaches the complementary burning line —0+ in the lower halfplane. 
There, a similar kind of deviation drives it to still another q(x"; t) to reach again 0+ and so 
forth and back. By subsequently running through these cycles between the caustics at ±0 + 
a net-motion into the direction of the barrier top may be generated. On the left side of the 
barrier (x, < 0) the same kind of scenario exists and at the top x = the burning lines 
intersect (depending in detail on V(q) within \q/l\ < 1, see fig [l].) so that small deviations 
in the vicinity of the bottleneck x = may lead from the set of right-barrier paths to that 
of left-barrier paths and vice versa. This allows for a motion starting in xi > to eventually 
reach the range on the opposite side of the barrier. So far the above discussion is restricted 
to class (c)-orbits with e rc = 0. However, for finite but small e re orbits merge close to the 
burning lines, and we find basically the same situation. The conclusion is that two real axis 
paths with TPs at q and — q , respectively, are linked by a sequence of real-time complex 
plane orbits tied together by small fluctuations near caustic lines. Since this under-barrier- 
motion is not a purely stationary one obeying (0), but can be seen as nearly stationary as 
it follows classical orbits most of the time, it describes quasi-stationary fluctuations (QSF). 
The QSF allow to move from > through the barrier range towards qj < and this way 
dominate in absence of true stationary points, 5S[q] = with q(Q) = q^ q{t) = qj, the path 
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integration in G(qf,qi,t) between the TPs. 

III. EXTENDED SEMICLASSICAL PROPAGATOR 

The action of the QSF can simply be approximated. For a cycle from Xi with q(xf,t) 
to x[ < Xi with q(x[]t) in the interval 5t we find with Cauchy's formula S(x' i , Xi, St) ~ 
— ESt where the short action is W(x',x) = J®, dq^j2m[E — V(q)] and the por- 
tion from the phase space deviation along 0+ is negligible. Accordingly, S(—qo, qo, At) « 
i\W(—q , q )\ — EAt where At is the time interval spent between the real axis TPs q , —go- 
Hence, one arrives at the crucial result that the real-time motion of the QSF gives rise to an 
imaginary part in the action which is identical to the known instanton or WKB exponent. 
The full action for a low energy motion from qi > to — qi now consists of two classical 
real axis segments from q { to q and from — q to —q i7 respectively, and QSF inbetween, i.e. 
S(—qi,qi,t) ~ 2W(qo,qi) + i\W(—qo,qo)\ — Et. In the semiclassical G(—qi,qi,t) the expo- 
nential of this action is accompanied by the contribution of Gaussian fluctuations around 
the real axis segments. 

The most powerful representation of the semiclassical propagator is the so-called 
Hermann-Kluk propagator (HK) || . It has the advantage of being determined by an initial 
value problem for the classical trajectories, namely, 

G HK (q f , qi ,t) = J^lh(q f , qi ,t,p,q)R(p,q,t)e lS ^ h (6) 

with the fluctuation prefactor R(p,q,t) and an overlap factor h(qf,qi,t,p,q) = 
(9/l7(P > 9 > *))(7(P,9)l9i) where 

(z| 7 (p, q,t)) = Q 1/4 exp{-| [x - q(t)] 2 + % -p(t)\x - <?(t)]} (7) 

is a Gaussian wave packet centered around the phase space point {p(t),q(i}}. In (|) one 
runs in the time interval t real trajectories from {p, q} to {p{t), q{t)} where the contribution 
of each orbit is weighted according to its action and fluctuation prefactor and the Gaussian 
overlap of its end-points with those of the propagator. 

The usual HK ([]) is exact for a pure parabolic barrier and therefore captures tunneling 
only for harmonic fluctuations around the barrier top but fails for long times || . To overcome 
this drawback we apply the results developed above and formally split the propagator in 
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phase space: G(qf,qi,t) = G > (qj,qi,t) + G<(g/, qi,t) where G> [G<] contains orbits with 
E >Vq [E <Vq — 5 p b] and G> coincides with Ghk ©■ Accordingly, G>(t) describes the 
time evolution up to moderate times (comprising the parabolic range Vq — <5 p b below the 
top) and G<(t) the long time behavior. Now, while in a strict sense the complex dynamics 
discussed above is only valid for very low E, we assume its applicability also for somewhat 
larger E and find with qi and qj on opposite sides of the barrier 

G < {q f ,q i ,t) = / ^h(q f ,q i ,t,P,q)R < ip,q,t)e iS< ^ t y n T(q ). (8) 

E<V ~8 pb 

Here, an orbit runs from {p, q} along the real axis to its TP {0, q }, jumps to {0, —go} to reach 
{p(t),q(t)} leading to a fluctuation prefactor R < and action S < . The position space jump 
costs T(q ) = exp[—\W(—q ,q )\/h} and eventually G< results from phase space averaging. 
In (H) we require \qi\/l, \qf\/l 3> 1 so that most of the dynamics is spent outside the barrier. 
Obviously, G < follows not just from switching in the integrand in G > to imaginary times in 
regions where E < V(q). 



IV. APPLICATIONS 



The extended HK (eHK) G e HK = C> + G < is now employed to scattering in an Eckart 
barrier Vb(q) = V / cosh 2 (g/Z) that has been of wide use, e.g. as a model for the H+H 2 
exchange reaction. Since asymptotically V(q) drops faster than any power of q we have 
<f>c — > and burning lines stretch parallel to the real axis. In fig. 0a the correlation function 
Cfi(t) = (ipf\exp(—iHt/h)\ipi) between two Gaussian wave packets is depicted. Initially, 
ipi [if)f\ is centered to the far right [far left] with V ^> pf/2 so that we are indeed in 
a deep tunneling regime. One clearly sees the exponential drop of G> and the startling 
accuracy of the eHK over the whole time range. The most sensitive observable for a real-time 
treatment is the transmission probability P(E) calculated by numerically Fourier transform 
Cfi(t). Remarkably, we get accurate data also for very low energies (fig. apart from small 



oscillations typical for real-time calculations [ lXf| . In the moderate energy range E/Vq > 0.5 
the "real-time" P(E) even improves the uniform WKB result. Here convergence for G< is 
achieved for roughly the same number of trajectories as in G> (typical number of trajectories 
for the set of parameters is 5 TO 4 ) so that in contrast to previous approaches [10] an extension 
of the eHK to two or three dimensional systems, which are of particular interest to study 
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chaotic tunneling, seems feasible. 

As an example where a dynamical approach is clearly needed we turn to an Eckart barrier 
driven by a periodic signal V = Vj, + qAsin(Qt) and focus on the range of non-resonant 
driving and weak to moderate driving amplitudes. In this case already the exact numerics 
is non-trivial since it is the long time tunneling behavior which is most sensitively affected 
by the driving and leads to a strong spreading of the wave packet. Typical results for the 
correlation function c/j(t) are shown in fig. §b. Compared to the static case one sees phase 
shifted oscillations and a revival type of phenomenon. Semiclassically, both effects originate 
from an intimate interference of (fast) above-barrier-paths (E > Vq), which cross the barrier 
and then are back-scattered, and (slow) driven tunneling orbits (E < Vq). Even in this 
time-dependent case the accuracy of the eHK is quite astonishing. 

V. CONCLUSION 

To conclude, our findings reveal for the first time how tunneling is encoded in the quantum 
propagator in terms of classical real-time orbits. In contrast to tunneling in the energy 
domain, real time barrier penetration cannot be described by individual tunneling orbits. 
Instead, it must be seen as a diffusion along a certain set of classical paths in the complex 
plane. This allows for a practical approach for semiclassical wave packet dynamics even in 
the deep tunneling regime of static and non-resonantly driven scattering processes. Explicit 
examples have been given for one dimensional cases, but the efficiency of the method suggests 
that at least two or three dimensional cases may be feasible. The main problem then will be 
that from a certain TP go a bunch of TPs on the other side of the barrier can be reached. 
This proliferation of orbits, however, seems tractable due to the exponential suppression 
of under barrier motion starting at go an d traveling over large distances. Work in this 
direction is in progress. Possible applications of our method and extensions of it in physics 
and chemistry may be e.g. mesoscopic systems in microwave fields or unimolecular reactive 
scattering. 
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FIG. 1: Orbits in the complex plane for V^g) = 1/(1 + q 2 ) 2 [thin; class (a) and (b) dotted, class 
(c) solid for various Xi and iji > 0, iji < 0]. Burning lines (thick) are shown for V^g) (dashed) and 
its asymptote l/q A (solid); dots are TPs. 

FIG. 2: Real part of eg vs. time for the static (a) and driven (b) scattering in an Eckart barrier. 
Parameters are 7Z 2 = 6, V /(p 2 /2) = 8, and (a) qi/l = —qj/l = 40, (b) qt/l = —qj/l = 15 with 
qi A/V = -0.75, n/^Vo/W = 0.02. 

FIG. 3: Transmission probability vs. E/Vq. Exact (solid), usual HK (dotted), eHK (dashed), and 
uniform WKB (dotted-dashed) are shown. 
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Figure 1 : Ankerhold et al 
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Figure 2: Ankerhold et al 




Figure 3 : Ankerhold et al 
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